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Abstract 

The free energies of six-vertex models on general domain D with various 
boundary conditions are investigated with the use of the n-equivalence relation 
which classifies the thermodynamic limit properties. It is derived that the 
free energy of the six-vertex model on the rectangle is unique in the limit 
(height, width) — > (00,00). It is derived that the free energies of the model 
on D are classified through the densities of left/down arrows on the boundary. 
Specifically the free energy is identical to that obtained by Lieb and Sutherland 
with the cyclic boundary condition when the densities are both equal to 1/2. 
This fact explains several results already obtained through the transfer matrix 
calculations. The relation to the domino tiling (or dimer, or matching) problems 
is also noted. 
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1 Introduction 



The six-vertex model is a solvable lattice model, usually introduced on the 
rectangle and considered with the use of the Bethe anzats method or the Yang- 
Baxter relation. The model was first solved by Lieb |ll3j and generally by 
S ut her land [1], both assuming the cyclic boundary condition in the horizontal 
and the vertical directions. In this paper /ls denotes the free energy obtained 
by Sutherland. 

We have, on the other hand, an examplej^j where the free energy is exactly 
obtained with another specific boundary condition. In this case the free energy 
is expressed in terms of the elementary functions, and thus it is apparently 
different from /ls- In the six- vertex model, the boundary condition is relevant 
even in the thermodynamic limit. This fact seems to be unusual comparing to 
other lattice models such as the Ising models. 

In this paper we consider the free energies of the six-vertex models in- 
troduced generally on domain D with continuous boundary and with various 
boundary conditions. The free energies are investigated and classified, with the 
use of the n-equivalence relation, and our study includes the cases where the 
transfer matrix method cannnot be directly applied. The main result of this 
paper is proposition 01 which states that the density of down arrows and that 
of left arrows on the boundary determine the free energy of the system on D. 
The free energy is identical to /ls when the two densities are both equal to 
1/2. This result also means that the free energy is still intensive even if the 
boundary effect remains relevant in the thermodynamic limit. 

Section T2.1 l is a short summary on the six-vertex model and the correspond- 
ing transfer matrix treatment. In section 12. 2 1 we introduce the domain D and, 
in section 12.3 I introduce an equivalence relation of boundary conditions called 
the n-equivalence 0. Two boundary conditions yield the identical free energy if 
they are n-equivalent. This n-equivalence is a generalization of the concept of 
boundary condition, classifies the infinite limit properties, and also corresponds 
to the irreducibility of the transfer matrix. In section 12.4 1 it is derived, with 
the use of the n-equivalence, that the free energy of the six-vertex model on the 
rectangle R with w columns and h rows is unique in the thermodynamic limit 
(w, h) — ► (oo, oo), specifically independent of the ratio w/h, independent of the 
order of two limits uu — > oo and h — > oo. The six- vertex model on a cylinder 
and a rectangle are considered and finally we obtain proposition 5. 

There exist several exact calculations which yield /ls with various boundary 
conditions. Our results explain why these free energies are equal to /ls and, 
in addition to it, can determine the exact free energies of six- vertex models 
which have not yet been solved. The results can also be written in terms of the 
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domino tiling language and also we can certify that proposition 5 is consistent 
with the results in this area. These facts concerning the relations with other 
known results are summerized in section 3. 

§ 2 The six-vertex models on D 
2.1 The six- vertex model 

Let us consider the square lattice and assign an arrow on each bond. The 
arrows are arranged such that two arrows come in and the other two go out 
at each site (the Ice rule). Then there exist six types of possible local arrow 
arrangements as shown in Fig.l. In this paper we are going to use the term 
'vertex' as a site and four bonds around it. Each vertex is assumed to have 
finite energy. The energy is assumed to be unchanged by reversing all arrows on 
the four bonds. Then we have three energy parameters and hence three types 
of Boltzmann weights a, b and c assigned to the vertices (see again Fig.l). We 
also introduce the Boltzmann constant fee, the temperature T, f3 = and 
the total number of sites N. The partition function is 

N 

z = e 

config. i=l 

where q is the energy of the i-th vertex and ^ config is taken over all the 
possible arrow configurations. The free energy / is obtained through —j3f = 
limiv^ooiV^MogZ. 

One can assign a line on each arrow pointing down or left (Fig.l). Then 
each arrow configuration corresponds to a line configuration on the lattice. The 
Ice rule corresponds to the restriction that each line begins from a bond on the 
boundary, continue until it reaches another bond on the boundary, and that 
the lines do not intersect each other. 

Let us consider a rectangle R with w columns and h rows. Assume that h 
and w are even. Let r\ = {x±, . . . , x m } be a line configuration with m lines on a 
row of vertical bonds in R, and let rf = {x[, . . . , x' m } be that on the row below. 
The symbol {x\,X2, ■ ■ ■ ,x m } denotes that there is a line on each x^-th bond 
(k = 1, . . . ,m) and otherwise there is not. The (?/, ^-element of the transfer 
matrix V is introduced as 

Vfj'rj = {x'i, . . . ,x' m \V\xi, . . . ,x m ) 

w 

config. k=l 
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where \x%, . . . ,x m ) is the state corresponding to rj = {x\, . . . , x m }, (k = 
1, . . . ,w) is the energy of the A:-th vertex on the row between the two rows 
of vertical bonds and here the sum X^config * s taken over all the possible line 
configurations with fixed r\ and rf . 

When we identify the bond at the right end with that at the left end on each 
row (we call this the cyclic boundary condition in the horizontal direction), the 
transfer matrix is identical for all rows. If the boundary condition is also cyclic 
in the vertical direction the partition function Z is written as 

Z = tvV h = ^ X i ~ \\ (h^oo), (2.3) 

i 

where Aj's are the eigenvalues of V and Ai > |Aj| for all i. 

Following the notations in [J] let us introduce A = (a 2 + b 2 — c 2 )/2ab. The 
transfer matrix is block diagonalized according to the number of lines m. The 
maximum eigenvalue of the transfer matrix lies in the block element with m = 
or w (A > 1), with m = w/2 (A < 1). When A > 1, we have two types of 
frozen phase where specific line configurations are dominant. In this case the 
free energy is a constant and all the arguments in this paper become trivial. 
We thus concentrate on the case with A < 1. 

2.2 Domain D 

We introduce the domain D. Let us consider a continuous and closed line 
7(t) = (x(t),y(t)) (0 < t < 1) which satisfies 7(^1) 7^ 7(^2) if t\ 7^ £2 except 
7(0) = 7(1). Let us assume that the sites are on the points (niai, n2<32), where 
n\ and 112 are integers, a\ and 02 are the lattice spacings. 

The sites inside 7 belong to D. The sites on the line 7 can be suitably 
defined to belong or not to belong to D. The vertices belong to D when the 
corresponding sites belong to D. Each row and column in D is assumed to be 
simply connected. 

The bonds/vertices in D are called the boundary bonds/vertices when they 
have non-zero intersection with 7. The sites of the boundary vertices are called 
the boundary sites. It is assumed that the number of vertices on the boundary 
divided by the total number of sites vanish in the thermodynamic limit. 

We will take the following limit: fix the line 7 and take the thermodynamic 
limit ai, (12 — * 0. This corresponds to taking the limit w, h — > 00 where w and 
h are the number of columns and rows in D. 
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2.3 The n-equivalence 



Now we introduce the n-equivalence [5]. Assume that each site and bond takes 
one of a finite number of states. Let us consider the set of sites which lie in D 
and can be reached from the boundary sites by n steps (n bonds) at minimum; 
we call these sites the n-boundary sites. Consider the set of bonds between 
(n — 1)- and n-boundary sites, and call them the n-boundary bonds. At last let 
us consider the n-boundary sites together with the n-boundary bonds and call 
them the n-boundary. Configurations on the n-boundary are called n-boundary 
configurations. Let {Fi} be the set of all the possible configurations on the n- 
boundary with a boundary condition V on the actually boundary of D. Two 
boundary conditions V and V are called n-equivalent when {r,} = {r^} as a 
set of n-boundary configurations. 

Proposition 1 Fix a sequence of finite lattices {D^}, where N is the number 
of sites and approaches to the thermodynamic limit as N — > oo. Suppose 
that the boundary conditions F and F' are n-equivalent, for each Dn, with 
n = o{N/N') where N' is the number of boundary sites. Then the two free 
energies with F and F 1 are identical in the thermodynamic limit. 

Proof: Let us write the partition function of the system on Dn as Z = ^ BiZ{. 
The factor Z{ is the partition function from the variables inside the n-boundary 
with a fixed n-boundary configuration Tj. The factor Bi is the contribution 
from the other variables with the boundary condition F and the fixed Tj. The 
factor BiZi is thus the partition function with the F and the IV From the 
assumption that n = o(N/N'), we obtain log Zi = —fiNfi + o(N) and log Bi = 
o(N). The index i runs from 1 to i max , where i max is the number of permitted 
configurations on the n-boundary and satisfies i mSLX < 0(r N ) where r is a 
constant. Then we obtain 





i=2 i=ii+l 



Ph (N — ► oo), 



(2.4) 



where /; = /i (1 < i < fi < fi (h + 1 < i < imax)- I 
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Note that what is important is the set of possible configurations, we do not 
need think about the number of ways in which each configuration is realized. 

This n-equivalence is a generalization of the concept of boundary condition, 
and can be considered in other lattice models such as the 19- vertex model, and 
also can be introduced in stochastic processes 0. All the boundary conditions 
are 1-equivalent, for example, in the case of Ising models with finite interac- 
tions, because the Ising spin states are independent of their nearest-neighbors. 
It follows that the free energies of the Ising models are independent of their 
boundary conditions. 

The equivalence can be introduced in a more generalized form. We can con- 
sider a boundary condition T on a subset of the boundary, and introduce the 
corresponding n-boundary and the corresponding n-equivalences. In proposi- 
tions |2 and |S] we concentrate on one and two of the four edges of the rectangle 
R and consider the corresponding n-equivalences. 

2.4 Results 

Lemma 1 Consider a rectangle R and assume the cyclic boundary condition 
in the horizontal direction. Then all the line configurations with m lines on 
the upper edge (the first row of vertical bonds) of R are 2m-equivalent to each 
other. 

Proof: Let {xi,X2, ■ ■ ■ ,x m }, Xi < a^+i, be a line configuration on the first row 
of R. Beginning from a line arrangement {xi, X2, £3, £4, X5} = {2,3,5,6,9} as 
shown in Fig. 2(a), for example, one can introduce the shift of lines {2, 3, 5, 6, 9} — 
{1, 2, 3, 5, 6} — ► {1, 2, 3, 4, 5}. When we generally have 

{1,2,..., i, Xj,Xj + \, . . . , x m _fc}, i+1 < Xj, after k-th step, one can introduce the 
line arrangement {1, 2, . . . , i, i + 1, xj, . . . , x m -(£:+i)} as the next one. We need 
m — I steps for the shift {x%, X2, ■ ■ ■ , x m } — > {1, 2, . . . , m} if X[ < m < and 
this means {x±, X2, ■ ■ ■ , x m } is at most m-equivalent to {1,2,..., m}. Then all 
the line configurations {x' 1; x' 2 , ■ ■ ■ , x' m } are again possible on the 2m-boundary. 
I 

Let us consider line configurations, on the left and right edges, with 1712 lines 
in the first h p bonds on the boundary, and rri2 lines in the next h p bonds, and 
so on. In this case one can introduce the line density p2 on the boundary as 
92 = m 2 /h p . 

Lemma 2 Suppose that the boundary configurations on the right and the left 
edges of R are identical and fixed with the density P2 = m 2 / h p where < P2 < 1 . 
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Then all the line configurations with m lines on the upper edge of R are m- 
equivalent where fh = am + a' (a, a' are constant). 

Proof: Let us introduce the shift of lines {1, 2, . . . , i, Xj, Xj+i, . . . , x m -k}— > 
{1,2, ... + l,Xj, ... , x m _( fc+1 )}, i + 1 < Xj, on each row without lines on 
the right and left boundary bonds, and {x±, x' 2 , ■ ■ ■ , x' m }— >{x' m , x' 1: . . . , x' m _ 1 } 
on each row with lines on the boundary bonds, as shown in Fig. 2(b). The 
line configuration will be arranged as {x±,X2, • • • , x m } — 2, . . . , m} within 
m! steps where m! = [m/ (h p — m2)]h p + h p . 

Next we consider an arbitrary line configuration {x'(, x 2 , ■ ■ ■ , x'^}, x'( < x'- +1 , 
and the shift {1,2,..., m}— >{x'{, x 2 , ■ ■ ■ , x'^}- Let us introduce {x'-y, x' 2 , ■ ■ ■ , x' m } 
— > {x[, x' 2 , . . . , x' m } on each row without lines on the right and left bound- 
ary bonds, {1, 2, . . . ,i, x" +1 , . . . , x'^ t } — > {x'^, 1, 2, . . . , i — 1, x'( , x'( + i, . . . , xj n _i}, 
as shown in Fig. 2(c), on each row with lines on the boundary bonds. Fi- 
nally the configuration {x'( } is possible on the m-boundary where 
fh < m' + m" and m" = [m/m2]h p + h p . | 

Boundary line configurations with which no configuration is admitted on 
the whole of the lattice should be excluded from our argument, because in this 
case the system cannot be a six-vertex model. We also assume the convergence 
of the free energy in the sequential limits h — > oo and w — > oo, and another 
sequential limits w — > oo and h — > oo. 

Proposition 2 Consider the six-vertex model on a rectangle R. Assume the 
cyclic boundary condition in the horizontal direction, or otherwise assume that 
the boundary configurations on the right and the left edges of R are identical, 
periodic with fixed period h p , with line density equal to p2 = rn,2/h p satisfying 
< p2 < 1. Assume that h is always a multiple of h p . Then, 

i) the transfer matrix of each row, or the product of the transfer matrices 
of sequential h p rows, respectively, is block- diagonalized according to the 
number of lines and each block element is irreducible, 

ii) the free energy is unique in the limit (w,h) — > (oo,oo), specifically the 
limit is independent of the order of two limiting procedures w — > oo and 
h — > oo, and also independent of the ratio w/h when one take h — ► oo 
with fixed w/h. 

Proof: First let us assume that the boundary condition is cyclic in the horizontal 
direction. The case with the fixed boundary condition can be treated similarly, 
and will be considered at the last of the proof. Let m be the number of lines on 
the upper edge (the first row of vertical bonds) of R. Then there are m lines 
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on every row of R because of the line conservation and the boundary condition. 
Thus the transfer matrix V is block-diagonalized according to m. Let V m be 
the block element of V with a fixed line number m. All the elements of V m 
are non-negative because they are sums of Boltzmann weights. From lemma ^ 
all the line configurations with m lines are 2m-equivalent. Hence there always 
exist allowed line arrangements on the lattice for any line configuration with 
m lines on the upper boundary (the first row of vertical bonds) and for any 
line configuration with m lines on the n-boundary with n > 2m + 2m, and 
thus we find that all the elements of V^ 71 are positive. (All the elements of 
V^ m are already positive, which is obvious from the proof of lemma Hence 
we find that V^ 71 is irreducible, because any matrix is irreducible if all of its 
elements are positive. It follows that V m is irreducible because V^ m cannot be 
irreducible if V m is not, this proves i). Then the Frobenius theorem works and 
we know the followings. There exists an non-degenerate eigenvalue \i(w) > 
such that Xi(w) > \Xi(w)\ where Xi(w) (i > 2) are the other eigenvalues of V m . 
We also know that there exists an eigenvector associated with X\{w) with all 
the elements being positive, i.e. the projections satisfy (x±, . . . , x m \max}> 
where |max) is the eigenstate associated with the maximum eigenvalue of V m . 
These results are valid for every finite m. 

The partition function Z with finite w and h is written as 

Z — (a?j , . . . , x m | V \xi,...,x m ) 

= c 1 X 1 (w) h + Y, CiXi ^ h ( 2 - 5 ) 

i>2 

where {x±, . . . , x m } and {x' x , . . . , x' m } are the line configurations on the upper 
and the lower edges of R, respectively, c\ = (x^, . . . , x^Jmax) (max|xi, . . . , x m ) 
is positive and independent of h and here we have assumed |max) is normalized. 
(The coefficients q are independent of h when V is diagonalizable. Otherwise 
V is expressed in Jordan form, Cj (i > 2) are asymptotically bounded by 
polynomials of h with finite degrees, and the following argument remains still 
valid.) We will consider the limit of 

-Pfh,w = ^~logZ = -\ogXi(w) + -rz'(h,w) (2.6) 
nw w h 

z\h,w) = -log Cl [l + J2-(^n (2-7) 

i>2 1 1 

The factor z'(h,w) is finite when h and w are finite, z'(h,w) depends on 
h only through the second term and so remains finite in the limit h — > oo 
with fixed w. When we take the limit h — > oo in ()2.6|l . the term (hw)^ 1 logZ 
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is convergent, because the term w -1 log Xi(w) is independent of h and the 
facter h~ 1 z'(h,w) converges to zero. Taking w — ► oo afterwards, one find that 
tu _1 logAi(w) is convergent because (hw)~ 1 log Z is assumed to be convergent 
in this limit. 

On the otherhand, when one first take the limit w — > oo in Q2.6|) . one find 
that the factor z'(h,w) is convergent, because (hw)" 1 log Z is assumed to be 
convergent and now it can be used that w~ l log \\(w) is also convergent in this 
limit. Therefore z'(h,w) is bounded by a factor C which is independent of 
w, because convergent series are always bounded. Again note that z'(h,w) is 
finite in the limit h — > oo, and therefore the factor C can be taken as a constant 
independent of w and h. 

Hence from (|2.6[) we have found that 

I - -- log Ai («;)|<£ (2.8) 
w h 

which means that the convergence of fh w in the limit h —* oo is uniform 
throughout w. This proves ii) when we consider the known fact that the uni- 
formity of the convergence yields uniqueness of the limit of double series. 

The last fact is well known but we are going to show a proof of it. We 
have shown that the convergence in the limit h —* oo is uniform: there exists a 
number foo,w which satisfies that, for every e > there is an integer ho(e) which 
is independent of w, such that \fh,w ~ f<x,v>\ < e f° r an h > ho(e). The free 
energy is convergent in the next limit w — > oo: there exists a number / which 
satisfies that, for every e > there is an integer wo(e) such that \ foo, w — f\ < e 
for all w > wo(e). Then for all h,w > max{/io(e), wo(e)} we have 

\fh, w -f\< \fh, w ~ foo,w\ + \foo,w ~f\< 2e, (2.9) 

which means 

t lim J h , w = f (2.10) 
(h,wj— »(oo,oo) 

as a double series. We assumed the convergence fh jW — > fh,oo ( w ~^ °°)) t nus 
taking w — > oo in (|2.9|) one obtains \fh,oo — f \ < 2e which means 

lim lim f h>w = f. (2.11) 

h— >oo ui— >oo 

Taking /i — ► oo also in (|2.9)1 one obtains \foo,w ~ f\ < 2e which means 

lim lim f h>w = f. (2.12) 

At last let us consider the case where the boundary line configurations 
on the right and the left edges of R are identical, fixed and the lines are lo- 
cated periodically with the period h p on the edges, with the line density being 



9 



P2 = mz/ hp satisfying < p2 < 1. Then one can introduce a transfer matrix 
V = V1V2 ■ ■ ■ Vh p where V\, ■ ■ ■ , Vh p are the transfer matrices of sequential h p 
rows of R, respectively. The partition function is expressed as a linear combi- 
nation of \i(w) h / hp , where Xi(w) are now the eigenvalues of V = V1V2 ■ ■ ■ Vh p - 
Then lemma El works and we obtain the same result. | 

Note that the 2m-equivalence of line configurations corresponds to the irre- 
ducibility of the block element V^ m , and hence to the irreducibility of V m . The 
matrix V is irreducible if the line configurations used as a bases for the matrix 
representation of V are n-equivalent to each other for some finite n. When the 
boundary is cyclic in the horizontal direction, Xi(w) is already known and we 
have lim^^oo w^ 1 log Xi(w) = — /3/ls- Hence the free energy itself is actually 
obtained in the proof of proposition [3 Here we show what is obtained from 
our formula using only the n-equivalence and the uniform convergence of the 
free energy. 

Lemma 3 Consider the six-vertex model on a rectangle R. Assume the cyclic 
boundary condition in the horizontal direction, and assume that the number of 
lines on the upper and the lower edges are m = m(w). Then all the boundary 
line configurations with the same m(w) yield the identical free energy in the 
limit (w, h) — > (00, 00). Specifically we have f = /ls when m = w/2. 

Proof: The line configuration on the upper edge is 2m-equivalent to arbitrary 
line configurations in the block element of V with m lines, and hence 2m- 
equivalent to the cyclic boundary with m lines. We take the limit h — > 00 
with fixed w, and obtain the free energy with (w,h) — > (w,oo). The resulted 
functions are the same for all of these fixed and the cyclic boundary conditions 
with m lines on the upper and the lower edges. Next taking the limit w — ► 
00, the free energies remain identical. In particular, we obtain the maximum 
eigenvalue of V and the known free energy fis when m = w/2. | 

When the system is cyclic in two directions, we know that the maximum 
eigenvalue of the the row to row transfer matrix lies in the block-element with 
w/2 lines in each row. Because of the symmetry of the system, the maximum 
eigenvalue of the the column to column transfer matrix also lies in the block 
element with h/2 lines in each column, otherwise we have contradictions. 

When the boundary is cyclic in the horizontal direction with fixed h/2 lines 
in each column, the 'alternate' line configuration with P2 = 1/2 with the period 
hp = 2 is possible on the right and left edges, hence lemma [2 works, and in 
lemma El /ls appears with the restriction that the number of lines is h/2 in 
each column. (Note that proposition ^ is valid when at least one necessary 
configuration is possible on the lattice.) 
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Lemma 4 In lemma /ls appears with the restriction that the number of 
lines is h/2 in each column. 

Let us consider a fixed boundary line configuration with a periodic pattern 
on the boundary of R. We assume that w and h are multiples of w p and 
hp, respectively and the line density is p\ = m\/w p on the upper and the 
lower edges, p2 = m^/Zip on the left and the right edges. The boundary line 
configuration on the upper edge is identical to that on the lower edge, and the 
configuration on the right edge is identical to that on the left edge. The limit 
will be taken with fixed w p and h p . In the proof of the next proposition we do 
not need the explicit form of \i(w). 

Lemma 5 With these conditions the line densities p\ and p2 determine the free 
energy of the six-vertex model on R: f = f(pi,p2)- Specifically /(1/2, 1/2) = 
/ls- 

Proof: First assume < p\ < 1 and < P2 < 1. Lemma |2] yields that all the 
line configurations on the upper edge with the density pi are n-equivalent with 
some n which depends only on w, and all the line configurations on the lower 
edge with the density p\ are also n-equivalent with the same n. The boundary 
effect is relevant to the limit only through p\ when we take h — > oo with fixed 
w. Next taking w — > oo we obtain the thermodynamic limit. This argument 
is also valid for the right and left edges with the density p2 taking the limit 
w — > oo with fixed h at first, and next h — > oo. The limit is unique because of 
the proposition |2j 

Consider the case with p\ = P2 = 1/2. The boundary configurations on 
the upper and the lower edges are n-equivalent to the cyclic boundary with 
w/2 lines with some n which depends only on w. Taking h — > oo with fixed 
w, we obtain the limit identical to that obtained with the cyclic boundary 
condition in the vertical direction with w/2 lines together with fixed boundary 
configurations with p2 = 1/2 on the right and left edges. Next taking w — > oo 
we obtain the thermodynamic limit, and from lemma 0] and proposition |2l the 
limit is identical to /ls- 

When p2 = the free energy is (1 — p\)£\ + pi£2, and when p\ = we have 
(1 — P2)ci + P2^2- As for the cases with the densities equal to 1, one can use 
the fact that f(pi,p2) = /(l — p\, 1 — p%) which comes from the symmetry of 
vertex energies. | 

Let us consider the sequence of vertical boundary bonds on the line 7. Let 
us assume a fixed boundary line configuration with a periodic pattern on this 
sequence with the line density p\ = rii/w p , i.e. n\ lines on every w p vertical 



11 



bonds on 7. Also assume a fixed and periodic boundary line configuration on 
the sequence of horizontal bonds on 7 with the line density p2 = n2/h p , i.e. 
n,2 lines on every h p horizontal bonds on 7. Assume that there exists a line 
on a horizontal bond on the left edge of D if and only if there is a line on the 
horizontal bond of the same row on the right edge of D. Assume also that there 
exists a line on a vertical bond on the lower edge of D if and only if there is 
a line on the vertical bond of the same column on the upper edge of D. The 
limit will be taken with fixed w p and h p . With these conditions we can derive 
the following: 

Proposition 3 The line densities p\ and P2 determine the free energy of the 
six-vertex model on D: f = f(pi,p2)- Specifically we obtain /(1/2, 1/2) = /ls- 

Proof: Let Dq be a rectangle of the width Ax and the height Ay sufficiently 
large satisfying D C D . The rectangle D is divided into small rectangles Ri 
of the width Ax and the height Ayi = Ah, where the lower edge of Ri coincide 
with the upper edge of Ri+\. Let D' = Uji^, where R[ is a rectangle of the 
height Ayi and of the width Axj, satisfying R[ C Ri and R[ C D, and satisfying 
that the number of columns in each R[ is a multiple of w p taking its maximum 
with the restriction that R[ C D. The sites in R[ and those on the edges of 
R'i, if any, are assumed to belong to R[, with the exception that the sites on 
R[ n R' i+ i belong to R\. The sum of the line numbers on the upper and the 
right edges of each R[ is generally equal to the sum of the line numbers on 
the lower and the left edges of R[, because of the line conservation property. 
Assume that the boundary line configuration on each edge of D' is identical to 
the corresponding boundary line configuration of D, i.e. there exists a line on 
a horizontal bond on a right edge of D' if and only if there is a line on the same 
row at the right edge of D, and so on for other edges. 

We have assumed in each R[ that the number of columns is a multiple of w p . 
In addition to it, we assume in each R[ that the number of rows is also a multiple 
of h p . The latter choice corresponds to introducing a sequence 02 — ► where 

each Q2 satisfies Ah = n'h p a2 with n' = 1, 2, It is sufficient to consider this 

specific sequence because the difference from the remaining cases of 02 vanish 
as iV — > 00 together with the ratio N' /N where N' is the number of boundary 
sites. 

Because of the identical line configuration on the right and left edges, the 
number of lines on the upper and the lower edges of R[ are the same. All 
the possible configurations on the upper edge are n-equivalent with some n 
independent of 02, all the configurations on the lower edge are also n-equivalent, 
and hence they yield the identical free energy in each R\ in the limit 02 — > 0. 
Next taking a\ — > 0, we find that each rectangle R[ yields /(pi, P2) which is the 
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free energy obtained in lemma [5J Thus we obtain fjyi = f(pi,p2) where fxj> is 
the free energy on D'. 

The result is independent of the ratio A/i/Axj and the ratio can be taken 
sufficiently small. The ratio of the energy contribution, the contribution from 
D\D' over that from D, goes to zero in the limit Ah — ► 0, because the line 
7, which determines the boundary of D, is continuous. Therefore we obtain 
|/d — /d'I < e for arbitrary positive e. | 

The result means that the free energy is still additive even in the situation 
where the boundary condition remains relevant in the thermodynamic limit. 

These results are applications of 6 to the case of the six-vertex model, 
and give sufficient conditions to have /ls- It should be noted that we didn't 
need diagonalize sequential transfer matrices but it was sufficient for us to 
consider the configurations which are n-equivalent, for the purpose to classify 
the thermodynamic limit properties. 

§ 3 Conclusion 

Our propositions explain several results already obtained and also able to de- 
termine the exact free energies of six-vertex models which have not been solved. 

One can introduce boundary conditions, such as the cyclic boundary con- 
dition, in which various boundary configurations are admitted. If our vertex 
energies and the temperature satisfy A < 1, and if we can find configurations 
being n-equivalent to those with p\ = P2 = 1/2, then the free energy is /ls- If 
we are in the parameter region with A > 1, and if one of the boundary config- 
urations with (pi,p2) = (0,0) or (1,1) is admitted, the system falls in trivial 
frozen phases. 

Owczarek and Baxter 8 solved (by the Bethe ansatz method) the six-vertex 
model with the cyclic boundary and a 'free' boundary condition, respectively, in 
two directions. Batchelor et solved (by the Yang-Baxter relation) the six- 
vertex model on the rectangle R with a specific boundary condition in which 
they have assumed that the arrow at one end of a row points right (left) if 
that on the other end of the same row points left (right), and also assumed 
the cyclic boundary condition in the vertical direction. One can find that the 
boundary configuration with the line density equal to 1/2 is realized with the 
restrictions in both [S] and jH], and hence it can be directly derived from our 
results that the free energies of these systems are /ls- Furthermore one can 
easily construct a bunch of n-equivalent cases which have not yet been solved 
and are extremely difficult to solve directly by the Bethe ansatz method or the 
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Yang-Baxter relation, but we now know that the free energies of all of these 
cases should be /ls- 

We would like to note that equivalences of boundary conditions in the six- 
vertex model is also investigated in ^Uj yielding /ls- 

One of the most interesting problems related to our results is the domino 
tiling (see jllll2j b The problem is to find the number of possible ways to cover 
a region completely using dominos (1x2 rectangles). It is of course equivalent 
to find the number of dimer coverings on a given lattice. 

Each configuration of domino is expressed in terms of a height function 
h(x). The correspondence is unique except an overall constant. More precisely 
we introduce h{x) as follows: color the squares in a checker-board pattern, and 
h{x) increases one in each unit moving anti-clockwise around black squares on 
the boundary of dominos, and decreases one around white squares, as shown 
in Fig. 3(a). Not all the regions can be tiled using dominos. The necessary and 
sufficient condition for tilability is written in terms of the height function as 
follows: \h(x) — h(y)\ < d(x,y) for all x and y on the boundary of the region 
where d(x, y) is the minimal number of steps moving from x to y with only 
black squares on its left. 

KasteleynflBj and independently Temperley and Fisher Jlj derived the num- 
ber of tilings on the m x n rectangle and obtained that in the thermody- 
namic limit it behaves exp(mnG/ir) where G denotes the Catalan's constant 

G = l/l 2 — 1/3 2 + 1/5 2 — 1/7 2 H . The number of tilings on the Aztec 

diamond (the 'square' rotated 7r/4 tiled by horizontal and vertical dominos) 
is also obtained|15j exactly as 2 71 ( ri+1 )/ 2 where n is the half of the diameter of 
the region. This is completely different from that of the m x n rectangle. We 
recognize that the number of possible ways of tiling strongly depends on the 
shape of the boundary. 

Cohn, Kenyon and Propp 16 showed a variational principle for the number 
of tilings: assume that the region is tilable and sufficiently 'fat', and assume 
that the slope of the height function (s,t) = (dh/dx, dh/dy) is asymptotically 
constant on the boundary, then the asymptotic number of tilings per domino 
is a function of (s,t). 

It is known [5ll7j that the number of possible domino tilings per domino is 
equal to the partition function of the six-vertex model with a = b = 1 and 
c = \^2, i.e. A = 0. The equivalence is obtained from the correspondence of 
configurations of dominos and vertices shown in Fig. 3(b). 

The boundary of the mxn rectangle has constant slope (0, 0) and, as a six- 
vertex model, line densities are p\ = 1 and p2 = on the upper edge. Mixing 
two boundaries shown in Fig. 3(c) periodically, we find that the line densities 
vary as p\ = 1 — e and p2 = + e while the slope remains (0, 0). Taking e = 1/2 



14 



we have p\ = p2 = 1/2. This modification can also be done for other three 
edges and we find, from proposition |3J that the number of tilings is obtained 
from the partition function by Lieb and Sutherland with a = b = 1, c = \/2 
that means A = 0. The free energy is (see for example [7]) 



Counting the residues at z = i, 3i, 5i, . . . on the imaginary axis, we find the 
Catalan's constant G and obtain that —/3f = 2G/n. The factor 2 corresponds 
to the fact that the number of vertices is equal to the number of dominos and 
twice the number of squares, and the result from our proposition |3] is consistent 
with that previously obtained by Kasteleyn and by Temper ley-Fisher. The 
limit m — > oo and n — > oo in the domino case is unique, which is also consistent 
with our proposition^ 

The six- vertex model with the domain wall boundary condition! 5. corre- 
sponds to the Aztec diamond. In this case the boundary line densities do not 
satisfy our condition, the free energy is not indentical to fis- 

Quite recently, S.Sheffield informed me his work|18j in which he derived a 
variational principle for the systems with gradient Gibbs potential. His argu- 
ment is based on the DLR condition, equivalent to the variational principle 16 
in the case of the domino tiling problems, and corresponds to proposition El in 
the case of the six-vertex model while the treatments of the thermodynamic 
limit is different. 

The author would like to thank Professor V. Korepin and T. Hattori for 
discussions. 
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Figure Captions: 



Fig.l: Six vertices, corresponding line configurations and their Boltzmann 
weights. 

Fig. 2: Shift of lines introduced in the proof of Lemma 1 and Lemma 2. 

Fig.3: (a) The height function for domino configurations, (b) Correspondence 
between vertices and domino tilings. The horizontal and the vertical lines are 
the lattice for the six- vertex model, while the lines rotated ±7r/4 are the edges of 
dominos. (c) Two boundaries with constant tilt (0,0) but different line densities. 
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